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Abstract 

The development of exotic new materials, such as metamaterials, 
has created strong interest within the electromagnetics (EM) commu- 
nity for possible new phenomenologies and device applications, with 
particular attention to periodicity-induced phenomena, such as pho- 
tonic bandgaps (PBG). Within this context, motivated by the fairly 
recent discovery in X-ray crystallography of "quasicrystals," whose 
diffraction patterns display unusual characteristics that are associ- 
ated with "aperiodic order," we have undertaken a systematic study 
of how these exotic effects manifest themselves in the radiation prop- 
erties of aperiodically-configured antenna arrays. The background for 
these studies, with promising example configurations, has been re- 
ported in a previous publication [V. Pierro et al, IEEE Trans. An- 
tennas Propagat., vol. 53, No. 2, pp. 635-644, Feb. 2005]. In the 
present paper, we pay attention to various configurations generated by 
Rudin-Shapiro (RS) sequences, which constitute one of the simplest 
conceivable examples of deterministic aperiodic geometries featuring 
quasirandom (dis)order. After presentation and review of relevant 
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background material, the radiation properties of one-dimensional RS- 
based antenna arrays are analyzed, followed by illustrative numerical 
parametric studies to validate the theoretical models. Design param- 
eters and potential practical applications are also given attention. 



1 Introduction 

In an ongoing plan of investigations PP, |2j, we have recently been concerned 
with the study of the radiation properties of structures with aperiodic order. 
Interest in this type of geometries is physically motivated by the recent dis- 
covery of "quasicrystals" [I], i.e., certain metallic alloys displaying X-ray 
diffraction properties that are known to be incompatible with spatial periodic- 
ity. In electromagnetic (EM) engineering, this has so far produced interesting 
applications in the field of photonic bandgap "quasicrystal" devices (see the 
brief summary and references in pQ). 

We have initiated our series of prototype studies with the aim of exploring 
the radiation properties of selected categories of aperiodic geometries, rep- 
resentative of the various hierachial types and degrees of aperiodic order in 
the "gray zone" that separates perfect periodicity from absolute randomness. 
To this end, we started in pQ with exploring the basic radiation properties of 
two-dimensional (2-D) antenna arrays based on certain representative cate- 
gories of "aperiodic tilings" (see [3], [BJ for an introduction to the subject). 
In parallel, we also began investigating the radiation properties of selected 
categories of 1-D aperiodic sequences, which are more easily amenable to 
explicit analytic parameterization and can thus provide better insight into 
the relevant wave-dynamical phenomenologies. Within this framework, we 
explored in |2| the radiation properties of 1-D quasiperiodic antenna arrays 
based on the Fibonacci sequence, which admit analytic parameterization in 
terms of generalized Poisson summation formulas. In this paper, we procede 
along the same route, but focus on the opposite end of the above-mentioned 
aperiodic-order "gray zone", namely, the quasirandom type of (dis)order. 
One of the simplest and most intriguing examples is provided by the so-called 
Rudin-Shapiro (RS) sequences [I| 0. These sequences were developed in a 
pure-mathematics context by H. S. Shapiro and W. Rudin during the 1950's, 
pertaining to some extremal problems in harmonic analysis. More recently, 
RS sequences have found important practical applications as signal process- 
ing tools for spread spectrum communication and encryption systems [TP] . 
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llj. Concerning their wave-dynamical properties, available studies are prin- 
cipally focused on the transmission characteristics of 1-D RS-based dielectric 
multilayers (see, e.g., [T2"]). 

In line with our stated interest in the radiation properties of aperiodically- 
parameterized antenna arrays, we concentrate in the present paper on 1-D 
RS-based antenna array configurations. Our implementation via analytic 
and numerical studies makes use of well-established spectral and statistical 
properties of RS sequences. 

Accordingly, the remainder of the paper is laid out as follows. Section 
12] contains the problem formulation, with description of the geometry and 
review of the basic properties of RS sequences. Section |H] deals with the 
radiation properties of representative classes of RS-based antenna arrays; in 
this connection, potential applications to omnidrectional and thinned arrays 
are proposed and discussed. Conclusions are provided in Section HJ 

2 Background and Problem Formulation 

2.1 Geometry 

The problem geometry, a 1-D array of N identical antenna elements with 
uniform inter-element spacing d, is sketched in Fig. ^ Without loss of 
generality, the array is assumed to be placed on the with radiators 

at x n = nd, n = 0, N — 1; in view of the rotational symmetry around the 
x-axis, we adopt the 2-D polar (r, 6) coordinate system. Neglecting inter- 
element coupling, and for an implied exp(jut) time dependence, the n-th 
antenna element complex excitation is denoted by 

I n = I n exp(-jk Q ndrj), (1) 

where I n is a real excitation amplitude which, by assumption, can take on 
only two values (denoted by black and white dots in Fig. [TJ) related to the RS 
sequence described below (Sec. 12.2(1 . Moreover, k Q = u^/eo^o = 27r/Ao is the 
free-space wavenumber (with Ao being the wavelength), and rj describes the 
inter-element phasing. For observation distances r ^> 2(Nd) 2 / Xq (Fraunhofer 
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region), the array is characterized by the "array factor" [T3] 

JV-l 

F s(0 = ^/nexp(j2vr<) 

n=0 
N—l 

= J2 In exp ti k ° nd ( sin 9 - v)} = F p( )> ( 2 ) 

n=0 

where £ = ^(sin^ — 77) / Ao is a nondimensional angular spectral variable. In 
what follows, we shall refer to |F S (£)| 2 as the "radiation spectrum", and to 
\F P (6)\ 2 as the "radiation pattern". The 1-D results discussed here can be 
extended to 2-D geometry-separable configurations via conventional pattern 
multiplication [T5j . 

2.2 Rudin-Shapiro Sequences 

Rudin-Shapiro (RS) sequences are two-symbol aperiodic sequences with quasir- 
andom character [7| 0. In what follows, their basic properties are briefly 
reviewed, with emphasis on those aspects that are directly relevant to the an- 
tenna array application in Sec. EH (see |I1]-[IE! ^ or more details). In its basic 
form, a RS sequence is generated from the two-symbol alphabet A a = {—1,1} 
via the simple recursive rule 

OC = 1, a 2n = «n, «2n+l = (-l)™«n- (3) 

Thus, for instance, the first ten symbols in the sequence are 

111-111-1111. (4) 

It can be shown that, in the limit of an infinite sequence, the two symbols have 
the same statistical frequency of occurrence [T7]. Throughout the paper, the 
sequence {«„} in flHJ) will be referred to as an alternate RS (ARS) sequence. 
This sequence can also be generated via a two-step procedure involving a 
substitution rule defined on a four-symbol alphabet A s = {A, B, C, D}, 

A^ AB, B -> AC, C -> DB, D -> DC, (5a) 

and a final projection II : A s — > A a , 

11(A) = 1, 11(B) = 1, n(C) = -1, 11(D) = -1. (5b) 
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Another possible construction algorithm is based on the use of the so-called 
RS polynomials [IT] , [""B] , [H"] defined by the recursive rules 

P m+ i(0 = PUO + exp (j'27r2 m O Q m (£), P = 1, (6a) 
Q m+1 (0 = P m (0-exp(j'27r2 m OQ m (0, Qo = 1. (6b) 

It can be verified that P m and Q m in (JHJ are (2 m — l)-th degree trigonometric 
polynomials in the variable exp(j27r£), 

2 m_ 1 

p m(0 = ^2?n exp (j27r<) , (7a) 

n=0 
2 m -l 

g n exp {j2im£) . (7b) 

n=0 

The link between RS polynomials and sequences stems from the polynomial 
coefficients p n and g n in ((7|). For P m -type polynomials, the coefficient se- 
quence coincides with the ARS sequence, i.e., p n = a n , n = 0, 2 m — 1. For 
Qm-type polynomials, the coefficient sequence (technically defined as ARS 
Golay-complementary [2*0] ) is obtained in practice by reversing the sign of 
the second half {n > 2 m_1 ) of the ARS sequence, 

n < n < 2™- 1 - 1, 

qn \ -a n , 2 rn - 1 < n < 2 m - 1. 1 ; 

The reader is referred to [TH], [21], [22] for more details and possible gener- 
alizations of ARS sequences. Antenna arrays based on ARS sequences are 
treated in Sec. 13.11 

Another RS sequence of interest for this investigation (see Sec. 13. 2|) is 
related to the ARS sequence in (jSJ) via the simple linear transformation, 

Pn = (9) 

It is readily verified that the sequence {(3 n } in (jUJ) is composed of symbols 
from the alphabet Ab = {0, 1}, and will accordingly be defined as binary RS 
(BRS) sequence. 

Despite the fully deterministic and deceptively simple character of ([2]) 
and Q, RS sequences exhibit some statistical and spectral properties that 
are remarkably different from those of typical deterministic sequences, and 
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similar in some respects to those of random sequences. In particular, the 
correlation properties of the ARS sequence in (JHj) resemble those of white- 
noise-like random sequences |18j . 

1 N ^ 

lim T? y] a n®n+u = &u0, (10) 
N^oo iv 4 
n=0 

where 5 pq is the Kronecker delta symbol, 8 pq = 1 for p = q, 5 pq = otherwise. 

From the spectral viewpoint, ARS sequences exhibit interesting extremal 
properties. It can be proved that their truncated Fourier series amplitude is 
bounded by [S] 

AT-l 



N < sup 

0<£<l 



^a n exp(j27r<) <(2 + y/2)y/N. (11) 

ra=0 

Note that while the lower bound in (fTTj) holds for any sequence {a n } gen- 
erated from the alphabet A a = {—1,1} [IB], the upper bound is peculiar of 
ARS sequences. For a general sequence {a n } generated from the alphabet 
A a , the upper bound in (fTTj) would be iV (with the equality sign holding, e.g., 
for periodic sequences) [TT], [TH]. On the other hand, for random sequences 
in A a , the upper bound typically behaves like y/ N log iV fTIj . 

From pijl . it can be proved that the Fourier spectrum of the ARS se- 
quence tends to a constant value in the infinite sequence limit [Hj, and is 
thus devoid of any strongly-localized spectral footprints. Again, this be- 
havior resembles the "flat-spectrum" character of white-noise-like random 
sequences. For the RS polynomials in © and (J7J), a behavior similar to (JTTj) 
is found, with a tighter upper bound [TTj 

|P m (0l<2^, |Q m (£)|<2^, (12) 

(compare with (fTTjl . with iV = 2 m ). The behavior of P m and Q m RS poly- 
nomials in a single period < £ < 1 is shown in Fig. for 1 < m < 4. It 
is observed that, at any order, the two polynomial types exhibit complemen- 
tary behaviors (maxima of P m correspond to minima of Q m , and viceversa). 
Moreover, all the odd-order (m = 21 + 1) RS polynomials have their only 
zeros at £ = 1/2 (P TO -type) and £ = 0, 1 (Q m -type) (see Figs. |2fa) and|^c)), 
whereas even-order (m = 21) polynomials do not vanish for < £ < 1 (see 
Figs. I2jb) andHfd)). Actually, it can be shown that, in general fTTj . 

PAO) = P 2 /(l/2) = 2\ P 2m (0) = 2 m , P 2m (l/2) = 0, (13a) 
Q 2 i(0) = -QM2) = 2\ Q 2m (0) = 0, g 2m (l/2)=2 m . (13b) 
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Concerning the BRS sequence in (j§J), its Fourier spectrum results from 
the superposition of a flat background and a periodic distribution (£ G Z) of 
spectral peaks (see [TJ], [23] for details). This could be expected by noting 
that the BRS sequence in Q can be interpreted as the superposition of 
an ARS sequence and a periodic sequence, whose cross-correlation can be 
proved to be zero [T7j . [THj . Quite remarkably, the BRS sequence exhibits 
the same Fourier spectrum as a Bernoulli random sequence in At = {0, 1} 
whose symbols are chosen with equal probability of occurrence [23], 
these two sequences are said to be homometric [2H] in technical jargon. 



3 Rudin-Shapiro Antenna Arrays 

We now investigate in detail the radiation properties of representative classes 
of antenna arrays based on RS sequences. Referring to the geometry in Fig. 
n the RS character is embedded in the array via the excitation amplitudes 
I n in (j2J). Two cases are considered: a) ARS sequence (J„ = a n ), and b) BRS 
sequence (/„ = (3 n ). Parametric studies have been performed, with emphasis 
on two meaningful array-oriented observables: the maximum directivity Dm 
and the side-lobe-ratio (SLR) ps, which can be compactly defined as 

2||^|| 2 , N 

Dm = -y^, (14) 



PS 



\F P f 

I PHoo 



In (fTlj) and (fT5|) . and ||-|| 2 denote the (supremum) and L 2 (rms) 

norms, respectively, 



«r/2 

\F p \\ 2 00 = sup \F P (9)\ 2 , \\F p \\l = / |F p (#)| 2 cosM?, (16) 



f7r/2 

~tt/2<8<tt/2 ~ J-tt/2 

and B denotes the main lobe angular region. Using one can show that 

IIP..2 o ,„V^ cos (kpvdri) sin (fc i/rf) 

||F P || 2 = 2c + 4 ^ c, — , (17a) 
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where the c v are correlation sums, 

N-l-v 

c v = J ^+- ( 17b ) 

1=0 

This provides the background for the theoretical considerations and syn- 
thetic numerical experiments that follow. 

3.1 Flat-Spectrum "Omnidirectional" Arrays 

We begin with the class of antenna arrays based on the ARS sequence {a n } 
in Q. A typical example of the radiation spectrum is shown in Fig. El for a 
100-element array. As could be expected from the spectral properties of the 
ARS sequence discussed in Sec. 12.21 the radiation spectrum exhibits a fairly 
flat behavior (no distinct localized spectral peaks) with truncation-induced 
rapid oscillations. This results in a radiation pattern with globally "omni- 
directional" characteristics. Interestingly, for this class of arrays, one can 
establish some general theoretical bounds on the maximum directivity Dm, 
which are practically independent of the inter-element spacing and phasing. 
First, one can show that the radiation pattern satisfies the constraint 

\\FpWl < (2 + V2) 2 N, \\F P \\ 2 2 « 2c = 2N. (18) 

In ()18)) . the upper bound on the norm follows directly from (fTTJ) . The 
approximation for the L 2 norm, numerically verified a posteriori (see below), 
is obtained from ()17aJ) by neglecting the ^-summation; this term vanishes 
exactly for d = p\o/2,p e N, and is more generally negligible since - in view 
of (HOI) and (fl7bjl . for arrays of finite but moderate-to-large size - one expects 
\cv\ ^ cq, for v > 1. Accordingly, it follows from (fTl| that the bound on the 
maximum directivity Dm is given by 

D M < (2 + V2) 2 « 10.7 dB. (19) 

By comparison with the computed value (via ([14))) in Fig. El (which is Dm ~ 
5.7, practically independent of d/\ and 77), this general bound turns out to 
be rather loose. In this connection, much tighter bounds can be established 
for arrays based on the RS polynomials in (jOJ). In view of (|7|1. these arrays 
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can be synthesized from the array factor F s (£) in (J2J) with J n = p n , g n and 
N = 2 m , yielding radiation patterns of the type 

IF /mi 2 - / \ p m[d(sm9 -ri)/\ }\ 2 , 

1 pWI "I |Q m [rf(sin^-r / )/A ]| 2 . (2Uj 

For this class of arrays, the relationships in (|18|) can be rewritten as 

\\F p \\l<2 m +\ ||F p ||^2-2 m , (21) 

where the upper bound on the norm now follows from (|12p. whereas the 
approximation for the L2 norm is obtained and justified as before (see ([18)1 . 
with AT = 2 m ). This yields the following bound on the maximum directivity 

D M < 2 « 3 dB, (22) 

which is consequently tighter than the general bound in (|19j) . By this syn- 
thesis, one could therefore obtain an antenna array with nearly "omnidirec- 
tional" characteristics. Moreover, for the special case of odd-order (m = 
21 + 1) RS polynomials, by exploiting the two available degrees of freedom 
(d/\ and 77) to scale and shift the spectral zeros in (|T3*j) . one can place up 
to two nulls at arbitrary directions, while keeping the global omnidirectional 
character. Two examples are illustrated in Fig. 01 in connection with 5-th- 
order RS polynomials (N = 32 antenna elements). Specifically, Fig. Efa) 
shows the radiation pattern synthesized with a P5 polynomial, with array 
parameters (d/X = 0.83 and 77 = 0.1) chosen so as to place two nulls at 
9 = —30° and 9 = 45°. Figure Efb) shows the radiation pattern synthesized 
with a Q5 polynomial, with array parameters (oI/Xq = 0.5 and rj — 0) chosen 
so as to place a null at broadside (9 = 0). The maximum directivity vs. 
the inter-element spacing is shown in Fig. |3Jc) for both configurations. An 
essentially flat behavior (magnified by the expanded scale chosen in the plot) 
is observed, with computed values very close to the theoretical bound in (|22|); 
this also serves as an indirect a posteriori validation for the approximation 
of the L2 norm used in (fTHjl and (j21j) . Clearly, in view of the omnidirectional 
character, with the absence of a main radiation beam, the SLR parameter 
ps in ()15|) is not of interest for this class of arrays. 

3.2 Thinned Arrays 

Next, we move on to antenna arrays based on the BRS sequence {f3 n } in (0). 
Assuming I n = (3 n in the array factor in (J2J), and recalling that (3 n G At = 
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{0, 1} (i.e., antenna elements are either "off" or "on"), one may interpret this 
choice as a special deterministic "array-thinning" strategy ^3] . Recalling the 
spectral properties of BRS sequences summarized at the end of Sec. 12.21 one 
can expect a radiation spectrum constituted of an ARS-type flat background 
(see Fig. EJ) with superposed discrete spectral peaks (grating lobes) typical 
of periodic arrays. These features are clearly visible in the example shown 
in Fig. E] for an array with N = 200 featuring N a = 90 "active" (i.e., 
I n = 1) elements. Based on these premises, and recalling the symbol equi- 
occurrence property for infinite RS sequences, one may anticipate radiation 
properties (in terms of beamwidth, directivity and SLR) comparable to those 
of a fully-populated periodic array, using nearly half of the elements, at the 
expense of a slight pattern degradation (due to the flat background spectrum) 
which is however expected (and acceptable) in most typical array-thinning 
strategies ^3] ■ in other words, for a fixed number of active antenna elements 
N a and array aperture La (defined as the distance between the first and 
last active elements), BRS-type arrays could yield an increase by a factor 
~ 2 in the maximum value allowable for the average inter-element spacing 
d av = L^/(N a — 1) so as to avoid the appearance, in the visible range, of 
grating lobes (d av ~ 2Ao vs. d av = Ao for periodic arrays). 

Figure |H1 illustrates the radiation characteristics of a typical BRS-type ar- 
ray with iV a = 90 active elements. Specifically, Figs. EK a) and Ef b) show the 
radiation patterns for average inter-element spacings d av = A and d av = 2A , 
respectively. Also shown, for comparison, is the reference behavior of a pe- 
riodic array with same number of elements and inter-element spacing (and, 
hence, comparable beamwidth). As expected, while the periodic array dis- 
plays the well-known grating lobes, the BRS-array does not. The expected 
slight degradation in the radiation pattern (due to the flat background spec- 
trum) is observed with amplitude below the secondary lobe level (~ — 13dB). 
This is better quantified in Figs. Etc) andEfd), where the maximum direc- 
tivity Dm and the SLR p$ vs. the average inter-element spacing are shown. 
It is observed that the BRS-array characteristics are only slightly worse than 
those of the periodic array for d av < Ao, and become almost uniformly better 
beyond the critical value d av = Ao at which the periodic array characteristics 
experience an abrupt degradation due to the entrance in the visible range 
of the first grating lobe (see Fig. Eta)). This trend holds up to values of 
d av ~ 2.2Ao, at which the first grating lobe of the BRS array enters the 
visible range. Qualitatively similar behaviors have been observed by us for 
arrays of different size. Some representative results are summarized in Ta- 
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ble I in terms of maximum directivity and SLR for the case d av = Xq and 
rj = 0. It is observed that, with the exception of small-size arrays (N a < 50), 
the SLR values are always comparable with the reference periodic-array level 
(~ -13dB). 

4 Conclusions and Perspectives 

This paper has dealt with the radiation properties of 1-D quasirandom an- 
tenna arrays based on Rudin-Shapiro (RS) sequences. Potential applications 
have been suggested, based on known spectral and statistical properties of 
these sequences. Results so far seem to indicate that RS sequences might 
provide new perspectives for the synthesis of antenna arrays. 

In particular, the ARS-type synthesis described in Sec. 13.11 might be of 
potential interest for smart reconfigurable antenna systems where omnidirec- 
tional receiving/transmitting features are needed, as a possible alternative 
to the code-based synthesis approaches in [23] • Other potentially interest- 
ing applications could be envisaged in the area of reflect-array synthesis of 
"simulated corrugated surfaces" j2S] and, more generally, "virtual shaping" 
applications in radar countermeasures j2Zj- m this connection, use could be 
made of recent advances in the practical synthesis of artificial impedance sur- 
faces j2H] for designing 1-D (strip) or 2-D (patch) planar arrangements based 
on suitable "dual" impedance surfaces. Examples are perfect-electric and 
perfect-magnetic conductors capable of providing the opposite-sign reflec- 
tion coefficients which, in a first-order physical optics approximation of the 
relevant scattering problem, play the role of the ARS excitation amplitudes 
in the array factor here. Such structures, currently under investigation, could 
yield strong reduction of specular reflection and uniform angular distribution 
of the backscattered power. 

On the other hand, the BRS-type synthesis described in Sec. 13.21 seems 
to offer a simple fully deterministic alternative to typical random thinning 
strategies for arrays of moderate-to-large size. In this connection, it is 
interesting to observe that, in view of the homometry properties mentioned 
at the end of Sec. El this strategy is in a sense equivalent to a completely 
random scheme, in which the antenna elements are turned on/off with equal 
probability. 

Planned future investigations involve the extension to 2-D geometries, as 
well as applications to reflect-array synthesis for radar cross-section reduction 
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and control. Exploration of the radiation properties of array configurations 
based on other categories of aperiodic sequences continues. 
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Table 1: BRS antenna array as in Fig. El Maximum directivity Dm and SLR 
Ps for d av = A , rj = 0, and various values of the number of active elements N a . 
The SLR values obtained from numerical simulation were almost independent 
of d av and rj. 



N a 


D M [dB] 


Ps [dB] 


10 


9.67 


-5.9 


25 


13.6 


-8.8 


50 


16.8 


-11.9 


100 


19.8 


-11.2 


250 


23.8 


-13.8 


500 


26.9 


-12.3 
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Figure 1: Problem geometry: A 1-D array of N identical antenna elements 
with uniform inter-element spacing d is considered, with antenna element 
positions at x n = nd, n — 0, N — 1, and uniform inter-element phasing rj. 
Black and white dots denote the two possible values of the real excitation 
amplitudes /„ in (JTJ), related to the RS sequence (see Sec. 12.2)1 . Also shown 
is the array-end-centered polar (r, 9) coordinate system. 
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Figure 2: Rudin-Shapiro polynomials P m (solid curves) and Q m (dashed 
curves), for various orders m. (a), (b), (c), (d): m = 1,2,3,4, respectively 
(see Sec. 12.21 for details). 



17 




Figure 3: ARS antenna array with N = 100 elements. Radiation spectrum 
|.F S (£)| 2 (normalized to its maximum value) vs. spectral parameter £ (see 
(J2J)). Maximum directivity is Dm ~ 5.7dB, almost independent of the inter- 
element spacing and phasing. 
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Figure 4: Examples of antenna arrays synthesized using the RS polynomials 
P5 and Q5 (N = 32 elements), (a): Radiation pattern \F P (6)\ 2 (normalized 
to its maximum value) for P 5 -type array, with d = 0.83Ao and 77 = 0.1 (with 
prescribed nulls at — —30° and 9 = 45°), having maximum directivity 
Dm = 2.9dB. (b): Radiation pattern for Qs-type array with d = 0.5Ao 
and 77 = (with prescribed null at 6 = 0), having maximum directivity 
Dm = 3dB. (c): Maximum directivity Dm vs. inter-element spacing for 
the two configurations in (a) and (b), shown as solid and dashed curves, 
respectively. 
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Figure 5: Thinned BRS antenna array with N = 200 (N a = 90 active ele- 
ments). Radiation spectrum |.F S (£)| 2 (normalized to its maximum value) vs. 
spectral parameter £. 
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Figure 6: As in Fig. 0E but radiation characteristics for t] = (solid curves), 
(a), (b): Radiation patterns |.F P (#)| 2 (normalized to their maximum value) for 
average inter-element spacings d av = Ao and d av = 2X , respectively, (c), (d): 
Maximum directivity Dm in (fT^j) and side-lobe ratio ps in (|T3j) . respectively, 
vs. average inter-element spacing. Also shown (dashed), for comparison, is 
the reference behavior of a periodic array of the same size. 
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